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Abstract 



We report results of a search for B ^ tv in a sample of 9.7 million charged 
B meson decays. The search uses both vrz^ and ivu decay modes of the r, 
and demands exclusive reconstruction of the companion B decay to suppress 
background. We set an upper limit on the branching fraction B{B — > rz^) < 
8.4 X lO"'' at 90% confidence level. With slight modification to the analysis 
we also estabhsh BiB"^ K^vi>) < 2.4 x 10~^ at 90% confidence level. 



1 



T. E. Browder/ Y. Li/ J. L. Rodriguez/ H. Yamamoto/ T. Bergfeld/ B. I. Eisenstein/ 

J. Ernst/ G. E. Gladding/ G. D. GoUin/ R. M. Hans/ E. Johnson/ I. Karliner/ 
M. A. Marsh/ M. Palmer/ C. Plager/ C. Sedlack/ M. Selen/ J. J. Thaler/ J. Wilhams/ 
K. W. Edwards/ R. Janicek/ P. M. Patel/ A. J. Sadoff/ R. Ammar/ A. Bean/ 
D. Besson/ R. Davis/ N. Kwak/ X. Zhao/ S. Anderson/ V. V. Frolov/ Y. Kubota/ 

S. J. Lee/ R. Mahapatra/ J. J. O'Neill/ R. Poling/ T. Riehle/ A. Smith/ 
C. J. Stepaniak/ J. Urheim/ S. Ahmed/ M. S. Alam/ S. B. Athar/ L. Jian/ L. Ling/ 
M. Saleem/ S. Timm/ F. Wappler/ A. Anastassov/ J. E. Duboscq/ E. Eckhart/ 
K. K. Gan/ C. Gwon/ T. Hart/ K. Honscheid/ D. Hufnagel/ H. Kagan/ R. Kass/ 
T. K. Pedlar/ H. Schwarthoff/ J. B. Thayer/ E. von Toerne/ M. M. Zoeller/ 
S. J. Richichi/° H. Severini/o P. Skubic/° A. Undrus/° S. Chen/^ J. Fast/^ 
J. W. Hinson/i J. Lee/^ D. H. Miller/^ E. L Shibata/^ I. P. J. Shipsey/^ V. Pavlunin/^ 
D. Cronin-Hennessy/^ A.L. Lyon/^ E. H. Thorndike/^ C. P. Jessop/'^ H. Marsiske/^ 
M. L. Perl/3 Savinov/^ D. Ugolini/^ X. Zhou/^ T. E. Coan/^ V. Fadeyev/^ 
Y. Maravin/^ I. Narsky/^ R. Stroynowski/^ J. Ye/^ T. Wlodek/^ M. Artuso/^ R. Ayad/^ 
C. Boulahouache/^ K. Bukin/^ E. Dambasuren/^ S. Karamov/^ G. Majumder/^ 
G. C. Moneti/5 R. Mountain/^ S. Schuh/^ T. Skwarnicki/^ S. Stone/^ G. Viehhauser/^ 
J.C. Wang/5 A. Wolf/5 j ^^15 g Kopp/^ A. H. Mahmood/^ S. E. Csorna/^ I. Danko/^ 
K. W. McLean/8 Sz. Marka/^ Z. Xu/^ R. Godang/^ K. Kinoshita/^'Q I. C. Lai/^ 
S. Schrenk/9 G. Bonvicini/o D. Cinabro/° S. McGee/° L. P. Perera/° G. J. Zhou/° 
E. Lipeles/^ S. P. Pappas/^ M. Schmidtler/^ A. Shapiro/^ W. M. Sun/^ 
A. J. Weinstein/i F. Wiirthwein/^'Q D. E. Jaffe/^ G. Masek/^ H. P. Paar/^ 
E. M. Potter/2 g Prell/2 y. Sharma/^ D. M. Asner/^ A. Eppich/^ T. S. Hill/^ 
R. J. Morrison/3 R. A. Briere/^ T. Ferguson/^ H. Vogel/^ B. H. Behrens/^ W. T. Ford/^ 

A. Gritsan/5 J. Roy/^ J. G. Smith/^ J. P. Alexander/^ R. Baker/^ C. Bebek/^ 

B. E. Berger/*^ K. Berkelman/'' F. Blanc/*^ V. Boisvert/^ D. G. Cassel/*^ M. Dickson/^ 
P. S. Drell/6 K. M. Ecklund/^ R. Ehrlich/^ A. D. Poland/*^ P. Gaidarev/^ L. Gibbons/^ 

B. Gittelman/6 S. W. Gray/^ D. L. Hartill/^ B. K. Heltsley/^ P. 1. Hopman/^ 

C. D. Jones/^ D. L. Kreinick/^ M. Lohner/^ A. Magerkurth/^ T. O. Meyer/^ 

N. B. Mistry/*' E. Nordberg/^ J. R. Patterson/*^ D. Peterson/^ D. Riley/^ J. G. Thayer/^ 
P. G. Thies/6 D. Urner/^ B. Valant-Spaight/^ A. Warburton/^ P. Avery/^ C. Prescott/^ 
A. I. Rubiera/^ J. Yelton/'^ J. Zheng/'^ G. Brandenburg/^ A. Ershov/* Y. S. Gao/® 

D. Y.-J. Kim/8 g^n^ Wilson^s 

^University of Hawaii at Manoa, Honolulu, Hawaii 96822 
^University of Illinois, Urbana-Champaign, Illinois 61801 
^Carleton University, Ottawa, Ontario, Canada KIS 5B6 

and the Institute of Particle Physics, Canada 
^McGill University, Montreal, Quebec, Canada H3A 2T8 

and the Institute of Particle Physics, Canada 



*Permanent address: University of Cincinnati, Cincinnati, OH 45221 

^Permanent address: Massachusetts Institute of Technology, Cambridge, MA 02139. 



2 



^Ithaca College, Ithaca, New York 14850 
^University of Kansas, Lawrence, Kansas 66045 
^University of Minnesota, Minneapolis, Minnesota 55455 
^State University of New York at Albany, Albany, New York 12222 
^Ohio State University, Columbus, Ohio 43210 
^''University of Oklahoma, Norman, Oklahoma 73019 
""^^ Purdue University, West Lafayette, Indiana 47907 
^^University of Rochester, Rochester, New York 14627 
^^Stanford Linear Accelerator Center, Stanford University, Stanford, California 94309 
^^Southern Methodist University, Dallas, Texas 75275 
Syracuse University, Syracuse, New York 13244 
^^University of Texas, Austin, TX 78712 
-"^^University of Texas - Pan American, Edinburg, TX 78539 
■"^^Vanderbilt University, Nashville, Tennessee 37235 
i^Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 
20 Wayne State University, Detroit. Michigan 48202 
^^California Institute of Technology, Pasadena, California 91125 
^^University of California, San Diego, La Jolla, California 92093 

^^University of California, Santa Barbara, California 93106 
^^Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
25University of Colorado, Boulder, Colorado 80309-0390 
2^ Cornell University, Ithaca, New York 14853 
^^University of Florida, Gainesville, Florida 32611 
^^Harvard University, Cambridge, Massachusetts 02138 



3 



The purely leptonic decay of the B meson offers a clean probe of the weak decay process. 
The branching fraction 



Stt V ^bJ 

exhibits simple dependence on the meson decay constant Jb and the magnitude of the 
quark mixing matrix element Vub- The dependence on lepton mass (m^) arises from helicity 
conservation and heavily suppresses the rate to light leptons. In the B system this means 
Tu is favored over fiu or eu final states. Nevertheless, the expected branching fraction 
B{B Tu) ~ 0.2 — 1 X 10~^ is small and the presence of additional neutrinos in the final 
state significantly weakens the experimental signature. 

In the context of the Standard Model, a crisp determination of CKM parameters may 
be obtained in principle by comparing B{B tv) with the difference in heavy and light 
neutral B^ masses [0, 

= ■^VBmBm'^f^BBSo{xt)\Vtd\^, 

a quantity which is known from Bd mixing measurements to considerable precision: 
Anid = 0.464 ± 0.18ps~^. In this comparison the dependence on the poorly known decay 
constant fs drops out, and one obtains 

B{B Tv) = ((4.08 ± 0.24) x 10"^) 

The range ±0.24 is set by current theoretical uncertainties. Given a sufficiently precise 
experimental measurement of the branching fraction, this relationship could be used to map 
out an allowed zone in the plane of Wolfenstein p and rj parameters that is roughly similar 
to that determined by measurements of \Vub\, but subject to a different mix of statistical, 
systematic, and theoretical uncertainties P]. Alternatively, if \Vub\ is obtained from other 
measurements in the B system, then the determination of B{B tv) may be viewed as a 
measurement of the decay constant Jb- This may be the only way to measure fs- Looking 
beyond the Standard Model, the B ^ tu rate is sensitive to effects from charged Higgs 
bosons and may be used to set a limit on charged Higgs mass. The sensitivity is greatest 
for large values of the Higgs doublet vacuum expectation value ratio, tan P . 

Existing experimental information is limited, however. A previous search by this col- 
laboration [0] in the T(4S') — > BB system yielded a 90% confidence level upper limit 
B{B Tu) < 22 X 10~^, and three searches |^ in the — > bb system have yielded 
upper limits ranging from 16 x 10~^ down to 5.7 x 10^^. Although the system offers 
powerful kinematical advantages, future measurements will be at the T(4S'). 

In this Letter we present results of a new search for B ^ tv using a method which 
is uniquely adapted to the T(4S') system. In this method we fully reconstruct the com- 
panion 5 in a quasi-inclusive reconstruction technique similar to that developed for earlier 
measurements 0. 

The data used in this analysis were collected with the CLEO II detector at the Cornell 
Electron Storage Ring (CESR). The data sample consists of 9.13 fb~^ taken at the T(4S), 



ub 



4 



corresponding to 9.66M BB pairs, and an additional 4.35 fb~^ taken below the BB threshold, 
which is used for background studies. 

CLEO II is a general purpose solenoidal magnet detector, described in detail else- 
where [^]. Cylindrical drift chambers in a 1.5T solenoidal magnetic field measure mo- 
mentum and specific ionization (dE/dx) of charged particles. Photons are detected using 
a 7800-crystal CsI(Tl) electromagnetic calorimeter covering 98% of An. Two-thirds of the 
data was taken in the CLEO II. V detector configuration, in which the innermost chamber 
was replaced by a 3-layer, double-sided silicon vertex detector, and the gas in the main drift 
chamber was changed from an argon-ethane to a helium-propane mixture. 

Track quality requirements are imposed on charged tracks, and pions and kaons are iden- 
tified by their specific ionization, dE/dx. Pairs of photons with an invariant mass within 2.5 
standard deviations of the nominal ir^ mass are kinematically fit with a 7r° mass constraint. 

mesons are identified in the Kg — ^ tt+tt" decay mode. Electrons are identified based 
on dE/dx and the ratio of the track momentum to the associated shower energy in the Csl 
calorimeter; muons over about 1 GeV/c momentum are identified by their penetration depth 
in the instrumented steel flux return; below about 1 GeV/c muons are not distinguished from 
pions. 

The experiment is fully simulated by a GEANT-based Monte Carlo |Tl| that includes 
beam-related debris by overlaying random trigger events on Monte Carlo-generated events. 
The simulation is used to study backgrounds and optimize selection criteria, but directly 
enters the analysis only through the calculation of the signal reconstruction efficiency. 

To search for B ^ ru decays we fully reconstruct each T(45') —>■ B^B~ event in the 
simultaneous decay modes B^ —> r^u ("signal -B") and B~ — > D^*^^{m:)~ ("companion 
5"). Here and throughout, charge conjugate modes are implied. 

For the signal B we accept any single track which passes track quality requirements. Pion 
candidates must have momentum greater than 0.7 GeV/c and must neither pass lepton iden- 
tification criteria nor be candidate daughters. We do not impose particle identification 
criterea. This approach encompasses the three decay modes r (e,/i)i/z/ and r — >■ vrz/, 
which together constitute 46.5% of the r branching fraction. Reconstruction efficiencies 
are 64%, 34%, and 84%, respectively, and there is some crossfeed into the "vri/" channel 
from the tau decay modes ez/z/, /iz/z/, and pz/. The crossfeed efficiencies are 6%, 20%, and 
8% respectively. The total r reconstruction efficiency, including r branching fractions and 
crossfeeds, is 32.9%. 

For the companion B, we take advantage of the large (46%) h cud branching fraction 
and seek to reconstruct B^ D^*^^(n7c)~, accepting either or D*'^ — i5°(7,7r°) and 
reconstructing the in the following eight modes, K~7t~^, K~tt'^tt^, K^ti^h^h^ , K^tt^tt^ , 
i^^7r"'"7r°7r°, K~tt~^tt^7i^tt'^, K^ti^tt^tt^ , and K'^n'^. Based on the reconstructed D° mass, 
the 7r° mass, and the kaon and pion particle identification information, we compute a 
quality factor and use it to reject poor candidates. The (nvr)" system may be any of the 
following: vr^, tt^tt^tt^, 71^71^71^71^11^ , 71^71^, 7r^7r^7r^7r'^, or tt^tt^tt^. 

With each B reconstructed in one of the target decay modes, we now require that there 
be no additional charged tracks in the detector, and that the sum of all energy in the 
crystal calorimeter not associated with the ionization energy deposition of charged tracks 
be less than a mode- dependent value Emax- For the clean decay modes of the companion B, 
B+ L)(*)07r+ and B+ D(*)°7r+7r°, we set E^ax = 0.6 GeV, while for all other modes it 



5 



is tightened to 0.4 GeV. The main source of non-associated calorimeter energy deposition is 
from hadronic interactions in the calorimeter that cast debris laterally and result in small 
energy deposits that are not matched with a parent track. Monte Carlo simulation and 
careful investigation of appropriate data samples indicates that on average such deposits 
sum to 240 MeV per (signal) event. Additional contributions arise from beam-related debris, 
averaging 26 MeV per event and concentrated in the far forward and backward portions of 
the calorimeter; and from real photons from incorrect signal reconstruction, which average 
10 MeV per event. In addition to this summed energy requirement, we also test whether 
any unassigned calorimeter signal can be paired with an already identified photon shower to 
form an object with invariant mass within 2.5 standard deviations of the vr" mass. If such a 
pairing can be made, the event is rejected. 

We suppress background from BB events by imposing requirements on the value of 
g^, the invariant mass squared of the nvr system. For most of the nir states we demand 

< 2.0 GeV^, but for the case mr = tt"*" no restriction is needed, and for nvr = 7t~^tt~7t~^71° 
we permit < 2.5 GeV^. 

Backgrounds arising from e~^e~ —>■ qq events ("continuum") are distinguished by a jetty 
topology. To suppress these backgrounds we compute the direction of the thrust axis of 
the companion B candidate and measure the angle 6 to the direction of the lepton or pion 
of the T candidate. For a signal event these directions should be uncorrelated and the 
I cos ^ I distribution uniform, while for continuum background the correlation is high and 
I cos^l peaks at 1. We require | cos6\ be less than 0.90 and 0.75 for r iuu and r — ttz/ 
candidates, respectively. Continuum background is more severe in the ttu mode and demands 
the tighter cut. Additional backgrounds from e~^e~ t^t~ are suppressed by requiring the 
Fox- Wolfram |13| moments ratio II2/H0 to be less than 0.5, which favors spherical topologies. 
Contributions from two-photon events {e~^e~ 'j^*^e~^e~) are negligible. 

The identification of acceptable candidates for the r daughter, the D^, and nvr system, 
together with the absence of extra tracks or significant extra neutral energy, marks the 
appearance of a signal candidate. We now characterize these candidates by the kinematic 
properties of the companion B, since there is no additional information in the lone r daughter 
track. In particular we use the total momentum Pb and energy Eb of the companion B, 
computed from momenta and energies of its daughter products. These raw quantities are 
then recast as the more useful beam-constrained mass M{B) = (-Eboam ~ PbY^"^ and energy 
difference AE = Eb — -E^beam variables. If more than one candidate is reconstructed in a 
given event, the one with the highest value of AE is selected. 

Figure |l| shows the distribution of events in the AE — M{B) plane for Monte Carlo BB 
background, Monte Carlo continuum background, Monte Carlo signal, and for the actual 
data set. The Monte Carlo background samples represent equivalent integrated luminosity 
of, respectively, three times and two times the actual data sample. The clustering of signal 
Monte Carlo events inside the M{B) signal region but around AE ~ —0.2 GeV is due 
to reconstructing i?^ D*^{n'K)^ as B^ D^{mT)^ . In such cases, the absence of the 
appropriate soft 7r° or 7 from D*^ decay lowers the candidate's total energy. Events in this 
satellite peak constitute 24% of the total signal yield. 

We select events whose M{B) falls within 2.5 standard deviations of the true B mass, 
and extract the signal yield by fitting the resulting AE distribution. The net signal efficiency 
including all secondary branching fractions for the analysis is £ = 0.69 x 10~^. The signal 
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fit shape is the sum of a narrow {a = 24 MeV) Gaussian centered at AE = for the 
primary signal yield, and a wide Gaussian {a = 115 MeV) centered at AE = —164 MeV 
for the D*^ satellite peak. The shapes and the relative normalization of these Gaussians 
are determined by Monte Carlo. Residual backgrounds are modelled by a linear distribution 
whose slope is determined by fitting the data lying outside the 2.5cr window in M{B). We fit 
the AE distribution by an extended unbinned maximum likelihood method to obtain the 
total yield of signal and background; the AE shape parameters are fixed by the procedure 
described above and are not varied in the fit. Figure 0a shows the final AE distribution 
of data inside the 2.5 standard deviation signal region of M{B)] six events remain after all 
selection criteria are applied. Figure |^ shows the fit shape with normalization as resulting 
from the likelihood fit; the central value of the fitted yield is 0.96 events. 

The background level is consistent with Monte Carlo expectations given the selection 
criteria and the size of the data sample. Figure ^ shows a comparison of the AE distribution 
for Monte Carlo events and data. To increase the yield for this plot we have released the 
restriction on leftover tracks, and here require exactly one extra charged track. These data 
events thus constitute in a sideband to the signal region. There are 71 such events in 
data, and 68 predicted by Monte Carlo. As evident in the figure, the Monte Carlo also 
reproduces the AE spectrum of these events very closely. Examination of Monte Carlo 
background events in the signal region itself shows (a) that the background is composed of 
approximately equal amounts of BB and continuum events; (b) that the background in the 
T ^ -Ku mode is dominated by continuum while the background in the r iuu mode is 
dominated by BB; and (c) about 75% of all background events, whether BB or continuum, 
have a K^^ present. Were it available, hadronic calorimetry would help suppress some of this 
remaining background. 

The branching ratio is related to the signal yield Ngig by B{B tu) = Nsig/N^^e where 
^BB = 9-66 X 10^ is the number of charged B mesons in the data sample and e is the efficiency 
as given above. We crosscheck the efficiency by conducting a separate analysis identical to 
this one in all key respects except that the tu target signal is replaced by D*^l~v whose 
branching fraction is large and well-measured. To ensure as much topological similarity to 
the TV case as possible, we restrict this ancillary analysis to the low-multiplicity sub-mode, 
D° ^ K-'n+. We find a yield of N{D*^l 

^)&Ata. = 43.1 ± 8.4 events in data, and compare 
this to the Monte Carlo result N{D*^i~i>)uc = 30.4±4.3 where the error is primarily due to 
uncertainties in the B~ — >■ D*^i~i> branching ratio 0. The discrepancy between these yields 
is 1.3(7. We adopt a conservative course, using the efficiency determined by Monte Carlo, 
and assigning to it a relative systematic error given by 6e/e = ^ (4. 3/30. 4)^ + (8.4/43.1)^ = 
24.1%. 

Figure Qd shows the likelihood function (£) plotted as — 2 In C/C^ax versus B{B — >■ 
Tu). Also shown is the result of convolving the likelihood function with the systematic 
uncertainty distribution of the efficiency (assumed to be Gaussian). The systematic error 
on the efficiency is dominated by the 24.1% discussed in the preceding paragraph, but 
also includes contributions from reconstruction efficiency uncertainty and uncertainty in the 
efficiency of the non-associated neutral energy cuts. In total, the relative systematic error 
on efficiency is 24.4%. We integrate the systematics-convolved likelihood function to obtain 
a 90% confidence upper limit Bgo on B{B -> tu) from 0.90 = J^'° C{B)dB/ Jq C{B)dB. We 
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find: 

B{B Tv) < 8.4 X 10" 



at 90% confidence level. This approach can be shown JTSl to be equivalent to the assumption 
of a flat Bayesian prior probability for B{B ru) and is known to yield a conservative 
upper limit. A frequentist approach based on generating Monte Carlo experiments gives 
B{B Tu) < 7.4 X 10"^ at 90% confidence level 

We also investigate the decay mode K'^vv There is currently no experi- 

mental information on this decay mode although limits on the related decays B — > X^z/z/ 
and B — >■ K*'^i>u exist ITH]. The search strategy is the same as described above, but we 



require that the lone track on the signal side fail lepton identification. The expected mo- 
mentum distribution of the charged track peaks at ~ 2.5 GeV/c, so we retain the 0.7 GeV/c 
momentum requirement previously applied to the pion candidate in the tcu mode. The re- 
sulting set of three K^uu signal candidates is a subset of the six tu candidates. They are 
marked by shading in Fig. |^. We perform the same unbinned likelihood fit as above and 
obtain a central value yield of 0.81 events. The efficiency of the K^uu is e = 1.8 x 10~^; 
we find B{B^ — >■ K^iyp) < 2.4 x 10~^ at 90% confidence level. The efficiency is calculated 
using the form factor model of Reference ITTI, but it changes only negligibly if instead we 



use 3-body phase space and a constant matrix element. We corroborate our result with 
an independent analysis, which is based only on counting events and yields an upper limit 
B{B^ K^w) < 6.6 X 10"^ at 90% confidence level. 

We have reported an analysis of 9.66 million charged B meson decays which results in a 
conservative upper limit on the branching fraction B{B tv) < 8.4 x 10~^. We also modify 
the analysis slightly to establish B{B^ K'^vv) < 2.4 x 10~^ The method used is optimized 
for conditions available at T(4S') experiments, and we anticipate useful application of the 
method to other rare decay modes with large missing energy. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada, the A. P. Sloan Foundation, the Swiss National 
Science Foundation, and Alexander von Humboldt Stiftung. 
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FIG. 1. Distributions in AE (vertical axis) and M{B) (horizontal axis), (a) BB Monte Carlo; 
(b) continuum Monte Carlo; (c) signal Monte Carlo; (d) data. The dashed lines delineate the signal 
region. 
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FIG. 2. Final results, (a) The six events fitted. Shaded entries correspond to candidates which 
are simultaneously K^uu candidates, (b) The fit shape with normalizations as resulting from the 
fit. (c) Distribution in AE of Monte Carlo (solid) and data (points) for events with exactly one 
extra charged track, (d) —2ln C/Cmax versus B{B —>■ tu). Solid: statistical errors only; dotted: 
systematic errors included as described in text. 
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